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Abstract 

We study a complex formation between the DNA and cationic 
amphiphilic molecules. As the amphiphile is added to the solution 
containing DNA, a cooperative binding of surfactants to the DNA 
molecules is found. This binding transition occurs at specific density 
of amphiphile, which is strongly dependent on the concentration of the 
salt and on the hydrophobicity of the surfactant molecules. We find 
that for amphiphiles which are sufficiently hydrophobic, a charge neu- 
tralization, or even charge inversion of the complex is possible. This is 
of particular importance in applications to gene therapy, for which the 
functional delivery of specific base sequence into living cells remains 
an outstanding problem. The charge inversion could, in principle, 
allow the DNA-surfactant complexes to approach negatively charged 
cell membranes permitting the transfection to take place. 



PACS.05.70.Ce - Thermodynamic functions and equations of state 
PACS.61.20.Qg - Structure of associated liquids: electrolytes, molten 
salts, etc. 
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PACS.61.25.Hq - Macromolecular and polymer solutions; polymer melts; 
swelling 
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1 Introduction 



In the last few years gene therapy has received significant attention both from 
the scientific community and from the general public. The development of 
new techniques for transferring genes into living cells allowed for the potential 



treatment of several diseases of genetic origin [fi|-||ll|. The central problem 
of gene therapy is the development of safe and efficient gene delivery system. 
Since both the DNA and the cell membranes are negatively charged, the 
naked polynucleotides are electrostatically prevented from entering inside 
cells. Furthermore, the unprotected DNA is rapidly degraded by nucleases 
present in plasma [IT]. 

Although, much effort has concentrated on viral transfection, non-viral 
methods have received increased attention. This is mostly due to the possible 
complications which can arise from recombinant viral structures, and the 
consequent risk of cancer. In the non-viral category, the DNA — liposome 
complexes have shown the most promise. Cationic liposomes can associate 
with the DNA segments, neutralizing or even inverting the electric charge 
of nucleotides, thus significantly increasing the efficiency of gene adsorption 
and transfection by cells. 

In this paper we present a model of DNA—amphiphile solutions. We find 
that in equilibrium, solution consists of complexes composed of DNA and 
associated counterions and amphiphiles. As more amphiphiles are added 
to solution, a cooperative binding transition is found. At the transition 
point a large fraction of the DNA's charge is neutralized by the condensed 
surfactants. If the density of surfactant is increased beyond this point, a 
charge inversion of the DNA becomes possible. The necessary density of 
amphiphile needed to reach the charge inversion is strongly dependent on 
the characteristic hydrophobicity of surfactant molecules. In particular, we 
find that for sufficiently hydrophobic amphiphiles, such as for example some 
cationic lipids, the charge inversion can happen at extremely low densities. 



2 The model 

Our system consists of an aqueous solution of DNA segments, cationic sur- 
factants, and monovalent salt. Water is modeled as a uniform medium of 
dielectric constant D. In an aqueous solution, the phosphate groups of the 
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DNA molecules become ionized resulting in a net negative charge. The salt 
is completely ionized, forming an equal number of cations and anions. Simi- 
larly the surfactant molecules are assumed to be fully dissociated producing 
negative anions and polymeric chains with cationic head groups. 

Following the usual nomenclature, we shall call the ionized DNA molecules 
the "polyions" , the positively charged ions the "counterions" , and the nega- 
tively charged anions the "coions" . To simplify the calculations, all the coun- 
terions and coins will be treated as identical, independent of the molecules 
from which they were derived. The DNA strands will be modeled as long 
rigid cylinders of length L and diameter a p , with the charge — Zq distributed 
uniformly, with separation b = L/Z, along the major axis. The cations and 
anions will be depicted as hard spheres of diameter a c and charge ±q. For 
simplicity we shall also suppose that each one of the s surfactant monomers 
is a rigid sphere of diameter a c with the "head" monomer carrying the charge 
+q. The interaction between the hydrophobic tails is short ranged and char- 
acterized by the hydrophobicity parameter x ( see Fig- !)• The density of 
DNA segments is p p = N p /V, the density of monovalent salt is p m = N m /V, 
and the density of amphiphile is p s = N s /V, where iVj is the number of 
molecules of specie i and V is the volume of the system. 

The strong electrostatic attraction between the polyions, counterions, and 
amphiphiles, leads to formation of complexes consisting of one polyion, n c 
counterions, and n s amphiphilic molecules. We shall assume that to each 
phosphate group of the DNA molecule can be associated at most one counte- 
rion or / < l max surfactants. This assumption seems to be quite reasonable in 
view of the fact that the electrostatic repulsion between the counterions will 
prevent more than one counterion from condensing onto a given monomer. 
On the other hand, the gain in hydrophobic energy resulting from the close 
packing of the surfactant molecules might be able to overcome the repulsive 
electrostatic interaction between the surfactant head groups, favoring con- 
densation of more than one surfactant on a given monomer (see Fig. 2). The 
/ amphiphilic molecules form a "ring" of radius a around the central negative 
monomer of the DNA (see Fig. 3). If we assume that most of the hydrocar- 
bon chain of the associated surfactants is hidden inside the DNA molecule, 
the maximum number of surfactants in a ring can be estimated from the 
excluded volume considerations, l max = 2na/a c , where a = (a p + a c )/2 is the 
radius of the exclusion cylinder around a polyion. 

At equilibrium, each site (monomer) of a polyion can be free or have 
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one counterion or a ring of I = l,...,l max surfactants associated to it. We 
define the surface coverage of counterions as p c = n c /Z, and the surface 
coverage of surfactant rings &spi = ni/Z, where n c is the number of condensed 
counterions and rt\ is the number of rings containing / surfactants. Each 
polyion has a distribution of rings containing from one to l max surfactants. 
We shall neglect the polydispersity in the size of the complexes, assuming 
that all the complexes have n c counterions and n s amphiphilic molecules - 
in rings of {pi} — with 

n s =J2 Zlpi • (1) 
1=1 

The total charge of each polyion is, therefore, renormalized from — Zq to 
—Z e ffq, with Z e ff = Z — n c — n s ||12||-||16||. ^From overall charge neutrality, 
the density of free cations is p + = p m + (Z — n c )p p , the density of free anions 
is p_ = p m + Ps, and the density of free surfactants is p{ = p s — n s p p . We 
shall restrict our attention to the limit of low surfactant densities, so as to 
prevent micellar formation in the bulk. 

The aim of the theory is to determine the characteristic values of n c , n s , 
and the surface coverage by rings {pi}- To accomplished this, the free energy 
of the DNA — surfactant solution will be constructed and minimized. 



3 The Helmholtz free energy 

The free energy is composed of three contributions, 

F ^complex ^electrostatic ^mixing • (2) 

The first term is the free energy needed to form the isolated complexes. The 
second term accounts for the electrostatic interaction between the counteri- 
ons, coions, surfactants and complexes. Finally, the third term is the result 
of entropic mixing of various species. 

To calculate the free energy required to construct an isolated complex 
composed of one polyion, n c condensed counterions, and n s condensed surfac- 
tants, we employ the following simplified model. Each monomer of a polyion 
can be free or occupied by a counterion, or by 1 < I < l max amphiphiles 
(see Fig. 2). Therefore, to each monomer i we associate occupation variables 
a c (i) and {o~i(i)}, which are nonzero if that particular monomer is occupied 
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by a condensed counterion or a ring with I surfactants, respectively. The free 
energy of N p isolated complexes can then be written as 



/Complex = -N p In E e"^, (3) 

where the sum is over all possible configurations of counterions and surfac- 
tants along a complex. For a particular configuration u, the energy can be 
expressed as the sum of three terms, E v = E% + Ei + E3 . The first one is the 
electrostatic contribution arising from the Coulombic interactions between 
all charged sites of a complex, 

q 2 * [-1 + ajs) + Ytr x + gc(j) + T!r x iMj)] m 

1 2^ D\r(i)-r(j)\ ' U 

where we have assumed that the only effect of association is the renormal- 
ization of the effective charge of each monomer. The second term E2, is due 
to hydrophobic interactions between the surfactant molecules, 

v z lmax (1 A- l'\ 

£2 = f E E ^^K(j), (5) 

where, in order to simulate the short-ranged nature of hydrophobic interac- 
tions, the first sum is constrained to run over the nearest neighbors. The 
hydrophobicity parameter \ is negative, representing the tendency of the 
two adjacent surfactant molecules to expel water. We can estimate its value 
from the experimental measurement of the energy necessary to remove an 



amphiphile from a monolayer and place it in the bulk [17 



The third contribution E3, accounts for the internal energy of each ring, 

Z Imax 

£ 3 = EE4)^ (6) 

i 1=2 

Ei is the interaction energy between I surfactants forming a ring. Each ring 
contains a maximum of l max sites, which can be occupied by surfactants. To 
each one of these sites we associate an occupation variable t(j'), which is zero 
if site j is unoccupied by a surfactant and is one if it is occupied ( see Fig. 
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3). The interaction energy of surfactants forming a ring can then be written 

as 

g2 Imax '\ ^ Imax 

E ' = W g 2os i 11 (,| i - + 2 g t(!)t0) • (7) 

The first term of Eq.(|7]) is due to electrostatic repulsion between the surfac- 
tant head groups, while the second is the result of attraction between the 
adjacent hydrocarbon tails. 

The exact solution of even this simpler sub-problem (i.e. evaluation of 
the sum in Eq.(|])) is very difficult due to the long ranged electrostatic in- 
teractions. We shall, therefore, resort to mean-field theory, which works 
particularly well for long-ranged potentials. Evaluating the upper bound for 
the free energy, given by the Gibbs-Bogoliubov inequality, and neglecting the 
end effects we obtain, 

P -^complex /3Np[fel ~\~ fhyd fring fmix] • (8) 

The first term, 

Imax 

Pfei = £S[-l+p c + Y,lPi} 2 -£SN p , (9) 

i=i 

is the electrostatic interaction between the sites along one rod and is related 



to Ex. S is expressed in terms of the digamma function 18 



S = Z[V(Z)-V(1)]-Z + 1, (10) 

and £ = (3q 2 / Db is the Manning parameter [|1^, pO[j . The second term in 
Eq. (|), 

Imax (j, _|_ rn \ 

(3f hyd = p x (z - 1) £ ^-J-p mPn , (11) 

n,m 

is the hydrophobic attraction between the rings inside a complex. The third 
term is the free energy due to the electrostatic and hydrophobic interactions 
between the surfactants forming a ring, 

O 1 7 _i_ i , Imax f-%"\/ I max 

pf rm9 = m ; +0 EV+^EV+ (12) 

1=2 ''max [—2 

Imax f 7 \ Imax / 7 \ / 7 \ 

£ Z Pl lln — + £ Z Pl l max 1 - — In 1 - — 

l = X \<>max ) i = x \ '"max J \ '•max J 
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where Uq ~ 0.25126591, and the reduced temperature is T* = ksTDa/ q 2 . 
Finally, the free energy of mixing for rings and counterions of a complex is, 



Pfn 



Z(l -p c - ^Pi) + ln(l -p c - J2pi) + %vMPc + 
i i=i 

i-max 

Z piliapi - Zp\ial max + (13) 
i=i 



Zpl r 



1--LW1 1 



where to be consistent with the expression (|T2"1), we have included a contri- 
bution to the free energy arising from the azimuthal motion of condensed 
counterions around the polyion, i.e. the last two terms of Eqn. ([l3|) . 

Once a cluster, constructed in isolation, is introduced into solution, it 
gains an additional solvation energy due to its interaction with other clusters, 
free counterions, free coions, and free surfactants. The electrostatic repulsion 
between the complexes is screened by the ionic atmosphere, producing an 
effective short ranged potential of DLVO form [El]]-p5[]. The electrostatic 
free energy due to interactions between various clusters can be estimated 
from the second virial coefficient, 



f3F c 



n r 



, 27rN 2 a 3 e- 2Ka 
VT*(Ka) 4 K 2 (Ka) 



(14) 



where (ko) 2 = inpl/T* and pi = a 3 [p p (Z — n s — n c ) + 2p m + 2p s ] is the 
reduced density of free ions. The free energy due to interaction between the 
complexes and free ions and surfactants can be obtained following the general 
methodology of the Debye-Huckel-Bjerrum theory fi~3|, 14}] , [ 26|1- [ 32| , 



(3F C 



N p (Z -n c -n s 



1 (a/L) 
T*(na) 2 



x 



-2 \n(KaK\ (ko) ) + I(kcl) — 



ina) 



with 



Una) 



xK 2 {x) 
K 2 {x) ' 



(15) 



(16) 
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where K n is the modified Bessel function of order n. The contribution to 
the total free energy arising from the interactions between the free ions and 
surfactants is given by the usual Debye-Huckel expression j36|,|27] 



/3F i 



V 



ln(l + na c ) — na c + 



,K(Xr 



(17) 



This term is very small and is included only for completeness. 

The last contribution to the total free energy Eq. (fj), results from the 
entropic mixing of the counterions, coions, surfactant and complexes, 



mixing 



~t~ -^m— ~t~ F s ~\~ F c . 



(18) 



The free energy of mixing is obtained following the general ideas introduced 
by Flory |33fl, 



/3F m+ 
(3F m . 
(3F S 

f3F c 



N m+ In <p m+ - N m+ , 
N m _ In m _ - N m _ , 
N s ln(0 s / n s ) — N s , 
' {Z + n c + n s ) 



iV p ln 



Z + n c + n s s 



-N p . 



(19) 



In the above expression m+ denotes free counterions, m— free coions, s free 
surfactant molecules, and c complexes. The 



7ip + fOc 

a 



6 

Tip* 



snp : 



6 

J* 
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a 

a 
1 



4(a/L) 



-{n c + n s s) 



(20) 



are the volume fractions occupied by the free counterions, coions, surfactants, 
and complexes, respectively. 
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4 Results and Conclusions 



The equilibrium configuration of the polyelectrolyte-surfactant solution is 
determined by the requirement that the Helmholtz free energy be minimum. 
Since F is the function of n s ,n c , and the surface coverage by rings {pi}, 
minimization of F implies that 

f)F BF lmax BF 

6F = ^-5n s + ^-5n c + £ °—5 Vl = . (21) 
on s dn c ^ dpi 

Using the constraint Eq. (|l|), Eq. (^T|) can be separated into Z ma z+1 equations, 

dF 

dn c 

and 



(22) 



OF OF 

—Zl + — = 0; l = l...l max . (23) 
on s dpi 

The system of equations (^) and ( U$) can, in principle, be solved numerically. 
However, for reasonable values of l max this requires a significant numerical 
effort. Instead of pursuing this brute force method, we note that to a rea- 
sonable accuracy, the surface coverage by rings, {pi}, can be approximated 
by an exponential distribution al . 

Pi = r r • 24 

We have checked that this is, indeed, a good approximation by numerically 
solving Eq. ( |2"3"D for an isolated complex. Using ansatz Q23|), the total free 
energy becomes a function of n c , n s , and a. For a fixed volume and number 
of particles the equilibrium corresponds to the minimum of Helmholtz free 
energy, 

dF 

^ = ' < 25 > 

dF 

^r ' (26) 

dF . 

fc = °- (27) 
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These are three coupled algebraic equations, which can be easily solved nu- 
merically to yield the characteristic number of condensed counterions, sur- 
factants, as well as the shape of the distribution of ring sizes (a). In Fig. 4 
we present a numerical solution of these equations. As a specific example we 
consider a cationic surfactant with an alkyl chain of s = 12 groups. In this 
case the hydrophobicity parameter can be estimated []30| to be in the range of 
X ~ —3, bksT . To explore the dependence of condensation on the hydropho- 
bicity of surfactant, we shall vary this value within reason. The density of 
monovalent salt and the DNA is taken to be 18 mM and 2 x 10 _3 mM, 
respectively. 

The resulting binding isotherms are illustrated in Fig. 4. The fraction of 
associated amphiphilic molecules f3 s = n s /Z, is plotted against the density 
of surfactant for a fixed amount of monovalent salt, p m . For small concen- 
trations of cationic surfactant, few amphiphilic molecules associate with the 
DNA segments. At the certain critical concentration, however, the system 
forms surfoplexes |3(J |32| - complexes in which the charge of the DNA 
is almost completely neutralized by the associated amphiphiles. If the den- 
sity is increased further, on average, more than one surfactant molecule will 
associate to each phosphate group, leading to charge inversion of the surfo- 
plexes. For highly hydrophobic surfactants the charge inversion can happen 
very close to the cooperative binding transition. We note that our theory 
predicts the binding transition to be discontinuous, this, most likely, is an 



artifact of the mean-field approximation |32 



We have presented a simple theory of DNA — surfactant solutions. Our 
results should be of direct interest to researchers working on design of im- 
proved gene delivery systems. In particular we find that addition of cationic 
surfactants leads to a strong cooperative binding transition. This transition 
happens far bellow the critical micell concentration. A further increase of am- 
phiphile density can result in the charge inversion of the DNA — surfactant 
complexes. This regime should be particularly useful in designing gene or 
oligonucleotide delivery systems. Until now most of non-viral gene-delivery 
systems were in the form of lipoplexes — complexes formed by DNA and 
cationic liposomes. To form the liposomes, however, is required a significant 
concentration of cationic lipid. Unfortunately, at high concentrations both 
lipids and surfactants are toxic to organism. Our model suggests that the 
charge inversion can be achieved with quite small concentration of cationic 
amphiphile, if it is sufficiently hydrophobic. This should reduce the risk of 
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unnecessary medical complications. 
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FIGURE CAPTION 



Figure 1. A cylindrical polyion of diameter a p , length L, and charge — Zq, 
surrounded by spherical ions of radius a c and amphiphilic molecules of s 
monomers. Each monomer of a macroion is free or has one counterion, or a 
ring made of I amphiphilic molecules associated with it. 

Figure 2. Schematic representation of a complex. Empty sites (monomers) 
(— ), sites with associated counterion (c), sites with / associated amphiphiles 
(si)- 

Figure 3. Ring composed of I surfactant molecules, l max = 15. 

Figure 4. Effective binding fraction of amphiphiles f3 s = n s /Z, as a function 
of amphiphile concentration p s . The concentrations of DNA and of added 
salt is 2 x 10~ 6 M and 18mM, respectively. The length of the DNA segments 
is 220 base pairs. The solvent is water at room temperature, so that £ = 4.17. 

Figure 5. Average size of rings in a complex (parameters the same as in Fig. 
4). 
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